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ABSTRACT 

n  This  memorandum  presents  a  method  of  fitting  equally-  spaced  data 
by  the  use  of  Chebyshev  polynomials.  The  method  is  explained  and 
tables  are  giv.en  to  facilitate  its  application.  The  tables  provide 
for  fitting  with  all  ordinates  adjusted,  fitting  with  *ero  error  at 
one  end  point,  cr  fitting  with  zero  error  at  both  end  points  as  the 
situation  may  demand. 
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INTRODUCTION 


The  advantages  of  using  Chebyshev  polynomials  in  the  polynomial 
approximation  of  a  function  are  well  known  (see  reference  1).  A 
description  of  the  simple  method  of  using  ordinary  differences  in 
the  fitting  of  equally- spaced  data  seems  to  have  been  omitted  by 
most  writers.  This  memorandum  makes  an  attempt  to  supply  such  a 
description. 


PROPERTIES  OF  CHEBY3HEV  POLYNOMIALS  USED  IN  FITTING  EQUALLY- SPACED  DATA 


The  Chebyshev  polynomial  Tn(u)  -  cos  (n  arc  cos  u)  has  the  well 
known  property  that  for  -1  <  u  <  1, 

-1  <  Tn(u)  <  l  , 

with  every  local  maximum  value  being  1  and  every  local  minimum  value 
being  -1. 

A  polynomial  of  degree  «,  Pn(x)f  can  be  approximated  in  the  inter¬ 
val  a  <  x  <  b  by  a  polynomial  of  degree  n  -  1  by  using 

-  (*)  +  »  (l) 

where  u  and  x  are  related  by  a  linear  transformation.  The  maximum 
magnitude  of  error  committed  by  using  Pfl_1  instead  of  Pn  will  be|^n| 
provided  «  lies  within  the  interval  -1  5  u  5  1  for  o  <  x  <  b. 

Differences  or  derivatives  can  he  used  in  finding  Kn  and  /*„_ x. 

The  advantage  of  using  differences  with  equally- spaced  data  is  that  the 
transformation  relating  u  and  x  enters  only  in  identifying  the  adjust¬ 
ments  made  to  the  ordinates  of  y  =  Pn(x )  to  give  pointB  which  lie  on 
the  curve  y  =  Pn^1(x).  Thus  the  needs  for  "shifted"  Chebyshev  poly¬ 
nomials  and  for  "normalizing"  the  argument  of  the  function  being 
studied  do  not  arise.  It  will  be  assumed  that  Newton  divided  differ¬ 
ence,  Lagrange,  or  graphical  methods  have  been  used  to  reduce  any 

unequally- spaced  data  and  that  the  data  are  therefore  available  for 
equal ly-spaced  arguments. 

With  data  available  for  x  =  a,  x  =  o  +  A,  x  =  o  +  ’2h,  .  .  . 

*  =  a  +  nh  =  6,  the  values  of  KnTni «)  can  be  computed  by  finding  the 
nth  difference  of  the  function  and  multiplying  by  the  correction 
factors  from  the  appropriate  table.  Subtracting  values  of  InTn(u) 
from  values  of  the  function  will  give  ordinates  for  a  polynomial 
curve  of  degree  n  -  1.  For  finding  a  lower  degree  approximation,  a 
set  of  *  equal ly-spaced  points  cn  y  -  can  be  obtained  by 

interpolation  and  the  process  repeated  using  the  s-point  table. 

The  number  of  steps  thus  taken  will  be  decided  from  the  degree  of 
the  polynomial  approximation  required  and  the  behavior  of  the 
successive  values  of 
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TABLES  OF  Tn  USED  IN  FITTING  WITH  ALT,  OPDINATES  ADJUSTED 


For  this  case  u  ~  -1  at  *  -  a  and  u  -  I  at  x  -  b. 


TABLE  1.  FIRST  DEGREE  FIT  (THREE  POINTS) 


=  y#  -  •2y1  +  y0  , 
f8(«)  =  '2w*  -  1  . 


tfT,  =  4  , 

Ka  =  0.'25  AaPa(x)  . 


TABLE  -2.  SECOND  DEGREE  FIT  (FOUR  POINTS) 


A'y  -  y9  -  3ya  +  3yt  -  y0  , 
f.  (u)  =  4ua  -  3«  . 


t 


3« 

'27  r,<M) 

-3 

->Z7 

-1 

•23 

1 

J23 

3 

'27 

9  \9T9  =  64  , 

I9  =  0. 140  625  A*P9(x) 
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TABLE  3.  THIRD  DEGREE  FIT  (FIVE  POINTS) 


A4y  =  y4  -  Ay$  +  6yt  -  4 yx  +  y0  , 
T4(u)  =  8«4  -  8uf  +  1 


>2u 

'2  T4(m) 

->2 

■2 

-1 

-1 

0 

2 

i 

-1 

•2  | 

•2 

TABLE  4.  FOURTH  DEGREE 

ABy  =  y8  -  5y4  +  10y,  - 
r„(«)  =  16a"  -  20m*  •? 


M 

f.(«) 

-1.0 

-1.00000 

-0.6 

0.07584 

-0.'2 

-0.84512 

0.2 

0.  84512 

0.6 

-0. 07584 

1.0 

1.00000 

=  12  , 

JT4  =  0.08333  .  .  .  A4F4(x)  . 


FIT  ( SIX  POINTS) 


10ya  +  5y*  -  y0  , 


A*r.  =  19.6608  , 

Kt  -  0.03086263  A* />.(*)  . 


♦ 
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try  -  y9  -  6ye  15y4  -  -20yg  +  15yt  -  6 yx  +  *0  , 


'2 


TABLE  6.  SIXTH  DEGREE  FIT  (EIGHT  POINTS) 

A7y 

-  y7  -  7ye  +  ’21y8 

-  35y4  +  35ya  -  *21  y2 

+  ?xi  -  y0 

T,M 

=  64a7  -  ll'2u6  + 

56u8  -  7u  . 

V 

823,  543  T7(u) 

7 

-823,  543 

5 

-539, '285 

3 

33, 933 

823,  543  A7T7  - 

41, '287,  680 

1 

694, 511 

*7  = 

0.01994646 

1 

-694,511 

3 

-  33,933 

5 

539, '285 

9 


(*) 


7 


823,  543 


TABLE  7.  SEVENTH  DEGREE  FIT  (NINE  POINTS) 


A9y  -  y8  -  8y7  +  ’28y9  -  56y8  +  70y4  -  56ys  +  '28 y2  -  8y,  +  y0  , 
re(u)  =  1'28«8  -  ’256u®  +  160u4  -  32ua  +  1 


u 

512  rfl(«) 

-1.00 

512 

-0.75 

449 

-0.50 

-'256 

-0.'25 

-’223 

0.00 

512 

0.'25 

-'223 

0.50 

-’256 

0.  75 

449 

1.00 

512 

51'2  AeTe  =  40, 320  , 

Kq  =  0.01269841  A  a/>8  (x) 


The  method  of  extending  these  tables  is  obvious.  Higher  degree 
Chebyshev  polynomials  can  be  obtained  from 

^n\  ~  Tn  (u )  _  * 
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TABLES  OF  Tn  TffSED  IN  FITTING  WIIH  ZgM  ERROR  FOR  x  -  a 


F\cc  this  case  u  s  for  s  -  a  and  w  =  i  for  x  -  bf  where  «_  is  the 
rpot  of  Tn(u)  =  O  which  ia  nearest  ti  -  -l.  Values  of  Tn («)  are  then 
commuted  for  equally- soeoed  argument p  in  the  interval  «_<«<!. 


i 


TABLE  >2 a.  SECOND  DEGREE  FIT  (FOUR  POINTS) 
Zero  error  for  x  -  a 


A3y  =  y8  -  3ya  +  3yt  -  y0  , 
f8(u)  =  4u8  -  3u  , 

«_  =  -  0.86602  54038  . 


u 

r.M 

-  0.86602  54038 

0.00000  00000 

-  0.  24401  69359 

0.67393  15713 

0.37799  15321 

-  0.91794  85072 

1.00000  00000 

1.00000  00000 

A*ra  =  5.77564  02355  , 

=  0. 17314  0978  A*/>8  (x)  . 
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TABLE  5a.  THIRD  DEGREE  FIT  (FIVE  POINTS) 
Zero  error  for  x  =  a 


A4y  =  y4  -  4y,  +  6yt  -  4yt  +  y0  , 
T4(u)  =  8«4  -  8uf  +  1  , 

«_  =  -  0.92387  95325  . 


V 

r4(«) 

-  0.92387  95325 

0.00000  00000 

-  0.44290  96494 

-  0. '26 149  35811 

0.03806  02338 

0.98842  81359 

0.  51903  01169 

-  0. 57456  05703 

1.00000  00000 

1.00000  00000 

A4T4  =  10.  27478  54210  , 

K4  =  0.09732  563349  A4/>4(x)  . 
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TABLE  4  a.  FOURTH  DEGREE  FIT  (SIX  POINTS) 
Zero  error  for  x  -  a 


A*y  =  yB  -  5y4  +  10yt  -  10y8  +  5yt  -  y0 
Tt(u)  =  16m8  -  >ao«*  +  5m  , 
u_  =  -  0.95105  65163 


M 

T.M 

-  0.95105  65163 

0.00000  00000 

-  0.56004  52130 

-  0. 16381  92620 

-  0. 17063  39097 

-  0.  75* 12  06893 

0. '21957  73934 

0.89431  87770 

0.60978  86967 

-  0.13694  54280 

1.00000  00000 

1.00000  00000 

A"rs  =  17.3700  -2549 

=  0.05757  04394  A8P§(x) 
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TABLE  5n  Firm  DEGREE  FIT  (SEVEN  POINTS) 

Zero  error  for  x  =  a 

A*y  =  y§  -  6yg  +  15 y4  -  ’20yB  +  15y#  -  6  yt  +  y0  , 
T#(«)  =  32»#  -  48«4  +  10m*  -  1  , 

*  =  -  0. 96592  58263  . 


1  u 

r.M 

1 

1  -  0.96592  58263 

| 

0.00000  00000 

1  -  0.63827  19219 

i 

0. 53024  02655 

t 

\  -  0.31061  72175 

0.31860  50447 

0.01703  70870 

-  0.99477  93213 

0.34469  13912 
• 

0.51470  51028 

0.67234  56956 

0.-28417  94180 

1.00000  00000 

1.00000  00000 

A *r#  =  '28.  5087  '2144  , 

JTt  =  0.03507  69825  AVb(x)  . 
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TABLE  6  a.  SIXTH  DEGREE  FIT  (EIGHT  POINTS) 


Zero  error  for  x  -  a 

ATy  =  y,  -  7ye  +  >21ys  -  35y4  +  35y,  -  '21ya  +  7yx  -  y0  , 
Tj(u)  =  64 u7  -  112u8  +  56«®  -  7 u  , 

«  =  -  0.97492  79100  . 


u 

t7  («) 

-  0.97492  79100 

0.0000  0000 

-  0.69279  53514 

-  0.7990  2858 

-  0.41066  '27928 

0.1783  6705 

-  0. 12853  02342 

0.7846  9717 

0.  15360  23242 

-  0.8817  1724 

0.43*73  48828 

0.0143  8894 

0.71786  74414 

0.6272  7220 

1. 00000  00000  , 

1.0000  0000 

=  45.896359  , 

K1  -  0.02178  7984  A7P,(*)  ■ 
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TABU??  OF  f„  OSSD  IK  FITmij  WITH  ZERO  ERROR  Kh'.  x  =  a  AND  FOR  x  =  b 


For  this  case  a  -  «_  for  x  -  a  end  «  =  for  *  =  6,  where  «_  end  u+ 
are  the  roots  of  r„v«)  =  0  nearest  «  ~  -1  and  «  =  1,  respectively. 
Values  of  r„(«)  are  then  oo^juted  for  equally- spaced  arguments  in  the 
interval  *_  <  «  <  *4  . 


t 


TABLE  3b 


THIRD  DEGREE  FIT  (FITE  POINTS) 


Zero  error  for  *  -  a  and  for  x  -  b 


A4y  =  y4  -  4ys  +  6yt  -  4yx  +  y0  , 
T4(u)  =  8u4  -  8us  +  1  , 

«_  =  -«♦=-  0. 92387  °5325  . 


u 

r4(«) 

-  0.92387  95325 

0.00000  00000 

-  0.46193  97662 

-  0.  34283  00857 

0.00000  00000 

1.00000  OOOOO 

0.46193  97662 

-  0.  34283  00857 

0.92307  95325 

0.00000  OOOOO 

A4  f4  =  8.  74264  06856  , 

K4  =  0.11438  19169  A *P4(x)  . 
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TABLE  4b.  FOURTH  DEGREE  FIT  (SIX  POINTS) 


Zero  error  for  *  -  o  and  for  x  -  b 


ABy  "  y§  -  5y4  +  10yf  -  10yt  +  5 yt  -  y0  , 
Ts(e)  =  16*B  -  '2Qu*  +  5u  , 

«  s  -  =  -  0.95105  65163  . 


tt 

T.M 

-  0.95105  65163 

0.00000  00000 

-  0.57063  39098 

-  0.  10501  08166 

-  0.19021  13033 

-  0.81740  21616 

0.19021  13033 

0.81740  '21616 

0.57063  39098 

0.  10501  08166 

0.95105  65163 

0.00000  00000 

A*  r,  =  15.  297  9351  , 

/,  =  0.06536  92995  A8/*,  (x)  . 
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TABLE  5b.  FIFTH  DEGREE  FIT  (SEVEN  POINTS) 

—  —  ■—■I— *—  -  -  -  -  -  -  -  -  -  .  -  - 

Zero  error  for  x  =  a  and  for  x  =  b 


A#y  =  yg  -  6 yB  +  15y4  -  SOy,  +  15yt  -  +  y0  , 

J*.(m)  =  32k*  -  48m4  +  18m*  -  1  , 

-  0.96592  58263  . 


u 

T.M 

-  0.96592  58263 

0.00000  ooooo 

-  C. 64395  05508 

0.49208  59669 

-  0.32197  52754 

0.38581  81249 

0.00000  00000 

-  1.00000  OOOOO 

0.32197  52754 

0.38581  81-249 

0.  64395  05508 

0.49208  59669 

0.96592  50263 

0. ooooo  ooooc 

A#ra  -  25  A695  1214  , 

=  0.03895  ^201  A*F#(*)  . 
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TABLE  fib.  SIXTH  DEGREE  FIT  (EIGHT  POINTS) 


Zero  error  for  x  -  a  ^nd  for  x  -  b 

A7y  =  y7  -  7y0  +  '21ys  -  35y4  +  35y3  -  -21y2  +  ly7  -  y0  , 
?'7(u)  =  64u7  -  ll‘2u6  +  56u8  -  7 u  , 

«_  =  -  «+  =  -  0.97492  79100 


u 

T7(u) 

-  0.97492  79100 

0.0000  0000 

-  0.69537  70786 

-  0.  7836  8357 

-  0.41582  62472 

0. 1391  8140 

-  0. 13627  54158 

0.8174  1285 

0. 13627  54158 

-  0.8174  1285 

0.41582  62472 

-  0. 1391  8140 

0.69537  70786 

0.7836  8357 

0.97492  79100 

0.0000  0000 

A7T7  --  40.-4017  1072  , 

X?  =  0.02475  1427  A 7P7(x)  . 
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TABLE'S  OF  CORRECTION  .FACTORS  USED  IN  FITTING 


Using  equation  (l), 

.  ^  »  (l) 

it  is  seen  tnat 

A*pn  =  kh  trrn  , 

from  which 

/„  =  (4"  *■„)'*  AV„  . 

Having  vo1ues  of  Tn  and  Anfn  from  the  preceding  tables, 

[(A*fn)-1  Tn(u)]  =  Cn 

can  be  computed.  This  correction  factor,  C„,  would  then  be  multiplied 
by  An?  to  give  the  amount  which  should  be  subtracted  from  ordinatas  on 
the  y  Pn(x)  curve  to  obtain  points  on  th©  y  ~  Pn_l (*)  curve. 


ft. 


TABLE  8. 


% 


CORRECTION  FACTORS  FOR  FIRST  DEGREE  FIT  (THREE  POINTS) 


X 

c> 

a 

0.-250 

a  +  h 

-  0.'250 

a  +  -2A 

0.'250 

xm 

0.-250 

-  A2PaC2(xj)  , 

Maximum  \P  2  -  Px  I  =  I Aa/>2  I  C2(xra)  • 
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Ij  rE  CORRECTION  FACTORS  FOR  SECOND  DEGREE  FIT  (FOUR  POINTS) 


- 

CL  _  _ 

O 

-  0. 1406  '2500 

0.0000  0000 

a  +  A 

0.1197  9167 

0.1166  8517 

a  +  2h 

-  0.1197  9167 

-  0.1589  3450 

a  +  sA 

0. 1406  '2500 

0. 1731  4098 

xm 

i 

0.  1406  '2500 

0. 1731  4098 

P7ixx)  =  P9(Xi)  -  A Vfc  {  }Xj  , 

l^sl  =  Maximum  iPg  -  P2  (  ~  lA8/*,  I  {  }  xm  * 

-  Cg  (all  ordinates  adjusted)  , 

<  ► 

=  CL  (zero  error  for  x  -  a)  . 

V  ✓  • 
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TABLE  10.  CORRECTION  FACTORS  FOR  THIRD  DEGREE  FIT  (FIVE  POINTS) 


0.0833  3333 


0.0000  0000 


0.0000  0000 


a  +  k 


-  0.0416  6667 


-  0.0254  5003 


0.0392  1356 


a  +  2A 


0.0833  3333 


0.0961  9939 


0. 1143  8192 


a  +  3* 
fl  +  Ah  - 


-  0.0416  6667 


0.0833  3333 


-  0.0559  1947 


0.0973  -2563 


-  0.0392  1356 


0.0000  0000 


0.0833  3333 


0.0973  2563 


0. 1143  8192 


M**)  =  {  )Xi  * 

If,  I  =  Maximum  I P4  -  P9 1  =  lA4P4l  {  }x  , 


C4  (all  ordinates  adjusted)  , 

C;  ( zero  error  for  x  -  a)  t 
C*  (zero  error  for  x  =  a  and  for  x  = 
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'AnT,F  11.  CORRECTION  FACTORS  FOR  FOURTH  DEGREE  FIv  (SIX  POINTS) 


x 


a 

a  +  h 
a  +  Zh 
a  +  3A 
a  +  Ah 
a  +  5A  = 


* 


c 

..rs 

C* 

-  0.  -508  6265 

0.0000  0000 

0.0000  0000 

. 

0.0038  5742 

-  0.0094  3115 

-  0.0068  6438 

-  0.0429  8503 

-  0.0435  3020 

-  0.0534  3219 

0.0429  8503 

0.0514  8632 

0.0534  3219 

-  0.0038  5742 

-  0.0078  8401 

0.0066  5438 

0.0508  6263 

0.057  5  7044 

0.0000  0000 

0.0508  6263 

0.0575  7044 

0,0653  6830 

P4(*t)  -  p.(xi)  -  AV,  {  , 

|X#I  -  Maximum  |P#  -  P4 1  -  |AB28 1  {  }x  , 


>  = 


j  - 


C9  (all  ordinates  adjusted)  , 

(aero  error  for  x  =  a)  , 

C8  (aero  error  for  x  =  a  and  for  x 


=  M  . 
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FACTORS  POT  FIFTH 


FIT  (SEVEN  POINTS) 


a 

0.0316  40625 

0.0000  00000 

0.0000  00000 

•  +  k 

0.0103  73264 

0.0185  99230 

0.0191  70055 

a  +  'Zk 

0.0142  79514 

0.0111  75704 

| 

0.0150  30209 

a  +  5* 

-  0.0316  40625 

-  0.0348  93860 

-  0.0389  56720 

a  4-  -4* 

0.0142  79514 

0.0180  54305 

0.0150  30209 

a  +  5* 

0.0103  73264 

0.0099  68157 

0.0191  70055 

a  +  6h  -  b 

0.0316  40625 

0.0350  76985 

0.0000  00000 

*■ 

0.0316  40625 

0.0350  76985 

0.0389  56720 

,(*j)  -  Pg(*j)  -  * 

|F#I  *  Maxima  IF,  -  Fgl  ®  lA*F#l  {  },  , 


C9  (all  ordinates  adjusted)  , 
C9  (aero  error  for  x  =  a)  , 

C9  (aero  error  for  x  -  a  and  x 


it 


TABLE  13.  CORRECTION  FACTORS  FOR  SIXTH  DEGREE  FIT  (EIGHT  POINTS) 


X 

i  c, 

C} 

CT 

a 

-  0.0199  46459 

0.0000  00000 

0.0000  00000 

I 

a  +  A 

-  0.0130  61645 

-  0.0174  09221 

-  0.0193  97287 

o  +  Zh 

0.0008  2186? 

0.0038  86258 

0.0034  44936 

a  +  3A 

0.0168  21265 

0.0170  96969 

0.0202  32135 

a  4-  4  A 

-  0.0168  21265 

-  0.0192  10841 

-  0.0202  32135 

a  +  5A 

-  0.0008  2186? 

0.0003  13506 

-  0.0034  44938 

a  +  6  A 

0.0130  61645 

0.0  136  66996 

0.0193  97287 

a  +  7A  =  b 

0.0199  4645S 

0.021?  8?983 

0.0000  00000 

0.0199  46459 

0.021?  8?983 

0.0247  51427 

/>.(*<)  =  M*,>  -  av,{  , 

If, I  -  Maximum  lP7  -  =  , 

r  ''j  -  Cr  (all  jrdinates  adjusted)  , 

y  -  Cf  (zero  error  for  x  -  j]  , 

J  -  C*  (zero  error  for  *  -  *  and  x  -  6). 
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NUMERICAL  EXAMPLF, 

In  order  to  illustrate  the  proposed  method,  Tables  11  and  12  will 
be  used  to  approximate 

y  =  in(l  +  *)  ,  (0  <  x  <  l) 

by  a  fourth  degree  polynomial.  Suppose  the  following  values  are  given. 


6x 

|  *"(l  +  x) 

0 

0.00000  0000 

1 

0. 15415  0680 

2 

0. -28768  -2073 

3 

0.40546  5108 

4 

0.51002  5624 

5 

0.60613  5804 

6 

0.69314  7181 

From  the  above  values 

A*i«(l  +  x)  -  -  0.00025  e428  . 

Turning  to  Table  12  and  remembering  that  a  =  0  and  6A  =  1, 

A*/>f  C,(a)  =  -  0.00000  81768  , 

A *P9  C,( a  ♦  h)  =  -  0.00000  26807  , 

A *P9  Cf(a  +  2A)  =  0.00000  36902,  etc. 

These  values  are  used  to  find  points  approximating  the  logarithmic 
curve  and  lying  exactly  on  the  curve  y  =  F#(x). 


6x 

M*) 

0 

0.00000  8177 

1 

0.  15415  3361 

2 

0.-28768  5763 

3 

0.4054  5  6931 

4 

0. 51002  9314 

5 

0.60613  8484 

6 

0.69315  5358 
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These  values  can  be  used  directly  to  obtain  a  fifth  degree  polynomial 
approximating  Ln( l  +  x).  To  obtain  a  fourth  degree  approximation, 
fifth  degree  interpolation  is  used  to  obtain  six  points  with  equally 
spaced  arguments. 


X 

0 

0.00000  8177 

0.  2 

0.18232  8444 

0.4 

0.33646  9893 

0.6 

0.47000  1120 

0.8 

1.58779  2299 

1.0 

0.69315  5358 

From  these  values 

-  0.00114  0176  . 

Using  Table  11,  a  -  0,  A  ~  0.2 

A bP6  C8(a)  =  -  0.00005  7992  , 

A 6Fb  C6 ( a  *  h)  -  0.00000  4390  , 

A V6  C6(a  +  Zh)  -  0.00004  9010,  etc. 

These  values  are  used  tc  fine  points  approximating  the  logarithmic 
curve  and  lying  exactly  on  the  curve 

y  =  M*) 


X 

P4(x) 

o 

0.00006 

6170 

0.2 

A  1  Q  ^  A 

A  —  V 

4046 

0.4 

0. 33651 

0903 

0.6 

0. 46995 

'7  '  -50 
•>  »  ,7 

0.8 

0.  50 7 7 9 

6697 

1.0 

0.693C9 

7365 

The  points  on  y  -  P4(x)  give  the  coefficients  in  P4  so  that 

ln( 1  +  x)  =  P4(x)  =  0.00006  6170 

+  0.99627  6293  * 

-  0.46644  2009  x* 

4  0.21862  2302  x* 

-  0.05542  5391  x4  , 

for 


0  <  x  <  1  . 


If  the  function  being  fitted  had  been  a  sixth  degree  polynomial,  the 
maximum  magnitude  of  error  would  not  exceed  IfJ  +  W',, !  -  0.00006  6170. 
It  can  be  expected  that  in  this  case  the  maximum  magnitude  of  error 
will  not  be  significantly  larger  than  this  amount. 


COMPARISON  OF  MAXIMUM  MAGNITUDES  OF  ERROR  FOR  DIFFERENT  TYPES  OF  FIT 

The  gain  in  accuracy  resulting  'Yoro  the  use  of  Chebyshev  poly¬ 
nomials  is  shown  in  Table  14.  The  «„  for  unadjusted  ordinates  is 
the  maximum  rsgnitude  of 


s ( s  -  l)(s  -  2)...  (s  -  11  i  l) 
1  •  2  •  3. .  • n 


in  the  interval  0  <  s  <  1. 

In  making  an  (a  -  l)st  degree  Chebyshev  fit,  the  interval  a  <  x  <  b 
is  divided  into  a  equal  parts  instead  of"  -  1  equal  parts  as  is  the 
case  for  a  point  interpolation  with  unadjusted  ordinates.  Therefore 


<  >-• 

=  Cn  (all  ordinates  adjusted) 

=  (zero  error  for  x  -  a) 

-  C"  (zero  error  for  x  =  a  and  x  =  b)  , 
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TABLE  14.  e„  =  *„  Any  =  MAXIMUM  MAGNITUDE  OF  ERROR  IN  n 
POINT  (xflt  xQ  +  h,  xQ  +  (n  ~l)A)  INTERPOLATION 


*n 


n 

Unadjusted 

Ordinates 

All  Ordinates 
Adjusted 

'2 

0, 1250 

0000 

0.0625  0000 

3 

0.0641 

5003 

0.0416 

6667 

4 

0.0416 

6527 

0.0263 

6719 

5 

0.0302 

6194 

0.0166 

6667 

6 

0.0234 

7346 

0.0105 

9638 

7 

0.0190 

1625 

0.0067 

8013 

8 

0.0158  8792 

0.0043 

6307 

y(a)  Unadjusted; 
Other  Ordinates 
Adjusted 

Kfl 

y(o)  and  y(b) 
Unadjusted; 
Other  Ordinates 
Adjusted 

0.0513  0103 

- 

0.0307  94439 

0.03619  11533 

0.0188  6468 

0.0214  1988 

0.0117  4721 

0.0130  4653 

0.0074  0610 

0.0084  1342 

op 


APPLICATIONS  OF  METHOD 


It  is  obvious  that  equation  (l)  can  be  used  in  many  situations  not 
discussed  in  this  memorandum.  Unequally- spaced  data  or  data  outside 
the  interval  a  <  x  <  b  could  be  used.  It  is  hoped  that  the  tables  of 
correction  factors  will  encourage  the  use  of  Chebyshev  fits  by  those 
who  have  previously  felt  that  the  extra  effort  outweighed  the  benefits. 
These  tables  cover  the  cases  most  frequently  encountered  in  making 
simple  approximations  at  the  desk.  For  machine  computation,  coding 
can  be  prepared  for  use  of  higher  polynomials  or  other  spacing  of 
arguments. 
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